Saturation and non-linear effects in difFractive processes 
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Through a direct implementation of the saturation regime resulting from the unitarity 
limit in the impact parameter representation, we explore various possibilities for the energy 
dependence of hadronic scattering. We show that it is possible to obtain a good description 
of the scattering amplitude from a hard pomeron provided one includes non-linear effects. 
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1 Introduction 

The most important results on the energy dependence of difFractive hadronic 
scattering were obtained from first principles (analyticity, unitarity and Lorentz 
invariance), which lead to specific analytic forms for the scattering amplitude as a 
function of its kinematical parameters — s, t, and u. Analytic S'-matrix theory 
relates the high-energy behaviour of hadronic scattering to the singularities of the 
scattering amplitude in the complex angular momentum plane. One of the impor- 
tant theorems is the Froissart -Martin bound ^ which states that the high-energy 
cross section for the scattering of hadrons is limited by 

where sq is a scale factor and /i the lightest hadron mass {i.e. the pion mass). As 
the coefficient in front of the logarithm is very large, "saturation of the Froissart- 
Martin bound" usually refers to an energy dependence of the total cross section 
rising as log^ s rather than to a total cross section equal to iQJ . 

Experimental data reveal that total cross sections grow with energy. This means 
that the leading contribution in the high-energy limit is given by the rightmost 
singularity in the complex-j plane, the pomeron, with intercept exceeding unity. 
In the framework of perturbative QCD, the intercept is also expected to exceed 
unity by an amount proportional to At leading-log, one obtains a leading 

singularity at J — 1 = 12 log 2 (as /tt). In this case, the Froissart-Martin bound is 
soon violated. 

But this is not the whole story, as there is another important consequence of 
unitarity, which is connected with the impact-parameter representation. 
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2 Unitarity bound 

Unitarity of the scattering matrix SS^ = 1 is connected with the properties of 
the scattering ampUtude in the impact parameter representation as it is equivalent 
at high energy to a decomposition in partial-wave amplitudes. As energy increases, 
the scattering amplitude in impact parameter can saturate the unitarity bound for 
some impact parameter bs- 

To satisfy the unitarity condition, there are different prescriptions. Two of them 
are based on particular solutions of the unitarity equation. 

For two-particle elastic scattering, the latter can be written 



(2) 



3(pi,P2,out|T|pi,p2,in) = ^ / d-i5\y^pr-^qr\\T^a\ 

7 •' \r=l r=l / 

The scattering amplitude in the impact parameter representation is defined as 

/•OC 

T{s,t)= / 6dWo(6A)/(6,s). (3) 
Jo 

with 

9/(6, s)<l. (4) 

and 

[3/(s,6)]'+[5R/(s,6)]'. (5) 
One of the possibilities is obtained in the JJ-matrix approach 0]: 

^(^'^)^ l-it/(.,b) - 
This solution leads to a nonstandard behaviour of the ratio ^ 

^ - 1, (7) 

ftot 

as s — > oo. For the highest energies reached so far {^/s = 2.0 TeV), such a ratio 
is 0.25. It is not too far from the standard value 1/2, but it is very far from the 
solution in the {/-matrix representation. 

The second possible solution of the unitarity condition corresponds to the eikonal 
representation 



r(s,t)= / 6d6Jo(&A)(l-exp(-x(s,6)). (8) 
Jo 

with t ~ A^. If one takes the eikonal phase in factorized form 

X{s,b)^h{s)f{b), (9) 

one usually supposes that, despite the fact that the energy dependence of h{s) can 
be a power 

/i(s)-s'^, (10) 
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the total cross section will satisfy the Froissart bound 

(Jtot < a log^ s . (11) 

We find in fact that the energy dependence of the imaginary part of the amplitude 
and hence of the total cross section depends on the form of i.e. on the s and 
t dependence of the slope of the elastic scattering amplitude. If f{b) decreases as 
a power of b, the Froissart-Martin bound will always be violated. In the case of 
other forms of the b dependence, a special analysis |S] is required. Note that the 
eikonal form does not correspond to a saturation of the amplitude: in this case, 
S/(s,6) reaches the black disk limit only asymptotically. Hence, the saturation of 
the black disk limit and the eikonal representation lead to different results for the 
scattering amplitude in the momentum transfer representation. 

Let us take, as an example, the hard plus soft pomeron model [S] [7j which 
includes two simple poles (a soft and a hard pomeron) to describe pp and pp scat- 
tering. In this case, the pp-elastic scattering amplitude is proportional to the 
hadron form factors and can be approximated at small t by: 

3T"(s,t) w k(s/so)'^e"'i"°s(''Ao) ^_;,2(s/so)'^e"^*i°s(«Ao)l ^2^^^ ^ (^2) 

where hi — 4.7 and /i2 = 0.005 are the coupling of the "soft" and "hard" pomerons, 
and ei = 0.0072, a'l — 0.25, and 62 = 0.45, a2 = 0.20 are the intercepts and the 
slopes of the two pomeron trajectories. The normalization sq will be dropped 
below and s contains implicitly the phase factor exp(— i7r/2). F^(t) is the square 
of the Dirac elastic form factor, which can be approximated by the sum of three 
exponentials jH! . We then obtain in the impact parameter representation a specific 
form for the profile function r(6, s) which we show in Fig. 1. One can see 
that at some energy and at small b, T{b, s) reaches the black disk limit. For one- 
pomeron exchanges, this will be in the region ^/s « 1.5 TeV. If one adds to the 
model 2-pomeron exchanges, the resulting r2 will saturate at y/s — 4.5 TeV. 

Saturation of the profile function will control the behaviour of (Xtot at super-high 
energies. Note that one cannot simply cut the profile function sharply as this would 
lead to a non-analytic amplitude, and to specific diffractive patterns in the total 
cross section and in the slope of the differential cross sections. Furthermore, we 
have to match at large impact parameter the behaviour of the unsaturated profile 
function. We have tried some specific matching patterns which softly interpolate 
between both regimes. The interpolating functions give unity in the large impact 
parameter region and force the profile function to approach the saturation scale bs 
as a Gaussian. 

If we take a single simple pole for the scattering amplitude, with an exponential 
form factor, the radius of saturation, and its dependence on energy, can be obtained 
analytically: 
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b GeV-i 



Fig. 1. The profile function of proton-proton scattering, (hard line and circles - at 
^/s = 2 TeV without and with saturation; dashed line and squares - at ^/s — 14 TeV 
without and with saturation; dash-dotted line - the eikonal form 1131 1 at ^/s — 14 TeV) 




Fig. 2. (Ttot calculated in different approaches. (The hard line - saturation regime; the 
long-dashed line - with the eikonal representation (1131 : the dash-dot-dotted line - with 
soft pomeron and unitarized hard pomeron ;7i; the dash-dotted line - with only a soft 

pomeron) 



where B is the average slope at small t. Hence the total cross section grows loga- 
rithmically at medium energies and grows like log^ s at very high energies. 

We show in Fig. 2 the possible behaviours of the total cross section at very high 
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energies, depending on the model and on the unitarization scheme. We also show 
there the result of a simple eikonalization, where we took 

Tppis, t = Q) = 2 j d?h[l~ exp(/icikG(s, b))] . (13) 

and where /leik = 1-2 was chosen so that the values trtot determined by (|13|l and 
by the saturation procedure are equal at ^/s = 50GeV, and G{s,b) is the Fourier 
transform of 3T°(s,t) given in (|12|l . 



3 Non-linear equations 

The problem of the implementation of unitarity via saturation is that the match- 
ing procedure seems arbitrary. Hence we considered a different approach to saturate 
the amplitude. It is connected with the non-linear saturation processes which have 
been considered in a perturbative QCD context ^] . Such processes lead to an 
infinite set of coupled evolution equations in energy for the correlation functions 
of multiple Wilson lines jll| . In the approximation where the correlation functions 
for more than two Wilson lines factorize, the problem reduces to the non-linear 
Balitsky-Kovchegov (BK) equation |lll I12| . 

It is unclear how to extend these results to the non-perturbative region, but 
one will probably obtain a similar equation. In fact we found simple differential 
equations that reproduce either the ?7-matrix or the eikonal representation. We 
can first consider saturation equations of the general form 

^^^.(A'), (14) 

with N the true (saturated) imaginary part of the amplitude. We shall impose the 
following conditions: 

(a) — > 1 as s — > oo, 

(b) dN/d^ ^ as s ^ c», 

(c) S{N) has a Taylor expansion in N, and considering the first term only gives 
the hard pomeron iVbarc — f{b)s'^. Similarly, we fix the integration constant 
by demanding that the first term of the expansion in s"^ reduces to A^baro- 

Inspired by the BK results, we shall use the evolution variable ^ — logs. If we want 
to fulfil condition (c), then we need to take S{N) = AN + 0{N'^). Conditions (a) 
and (b) then give S{N) = A{N — N'^) as the simplest saturating function. The 
resulting equation 

^^A{N-N^) (15) 
a logs 

has the solution 

N= (16) 
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One can in fact go further: eq. p5|l has been written for the imaginary part of the 
amphtude. If we want to generahze it to a complex ampUtude, so that it reduces 
to (|15|) when the real part vanishes, we must take: 

8 A 

A{A + iA^). (17) 



(9 logs 



The solution of this is exactly the form (j^ obtained in the [/-matrix formalism, 
for 5[/(s,6) = s^f{b). 

Many other unitarization schemes are possible, depending on the function J-{N). 
We shall simply indicate here that the eikonal scheme can be obtained as follows: 

dN 

^ = A(l-iV)(-log(l-7V)). (18) 

Other unitarization equations can be easily obtained via another first-order 
equation. The idea here is that the saturation variable is the imaginary part of the 
bare amplitude. One can then write 

dN „ dN 9iVbarc 



= :f'{n)^- — = ^^^'(iv), (19) 

OA^barc 9 logs (9 log S 

with A^barc the unsaturated amplitude. 

This will trivially obey the conditions (a) - (c) above, and saturate at iV = 1. 
Choosing T'{N) ~ 1 — N gives the eikonal solution 

N{b, s) = 1 - exp(-iVbarc(6, s)) , (20) 

whereas J^'{N) = (1 — N^ leads to the [/-matrix representation 

We can come to the same results if we solve this equation via an iteration 
procedure. For that let us take some model of the hadron interaction in which the 
main hadron-hadron interaction is created by the valence quarks surrounded by 
clouds of sea quarks. 

Ref. pSI proposed a picture in which the whole impact parameter interval is 
divided into small and large distances: the BFKL approximation works within a 
small radius i?o and the BK representation works for large impact parameters. 

We examine two cases: the first includes only a short-range form factor (as a 
simple Gaussian), and the second uses a short-range form factor and a long-range 
one (in the form of a MacDonald function) . 

So, we divide the whole energy interval in small pieces inversely proportional to 
s^^. We then obtain the step in s 

so = s^/s^^ = 1/s^ (21) 

and the number of such pieces will be 

n = int {s^^) , (22) 
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where the function int chooses the nearest natural number. 
N{s,b) after this small energy interval will be 

No{so,b) = sof{b) = f{b)/s^. (23) 

We calculate the derivative 6 on this small interval and increase N{so, b) by 5Nq{sq, b): 

(24) 



5i = l-Ni, 



Ni+i = Ni + SNo{so,b). 

We then iterate the above procedure for Ni = N{si, b) until wc get to the end of the 
interval. It is clearly understood that, if the energy is sufficiently high, we obtain 
for some iteration Ni- = 1 and 5 = 0, so we reach the saturation bound at some 
impact parameter. Of course this energy will depend on the form of f{b). 
In this case 

T{s,t^ 0) = / bdbN{s,b) (25) 
Jo 

and 

<Ttot(s)=47rS>r(s,t==0). (26) 

The results of our calculations are shown in Fig. 3. We can see that energy depen- 
dence of the total cross section is a simple logarithm. So, in this case, (Ttot("S) docs 
not saturate the Froissart unitarity bound and coincides with the eikonal solution 
with a Gaussian eikonal. 

Note that wc obtain in this case the right energy dependence not only for (Ttot(s), 
but also for the differential cross sections. Of course, this picture is only qualitative. 




Fig. 3. (Jtot calculated in different approaches. (The hard line - f{b) in Gaussian form; 
the long-dashed line - with the contribution of the large distances; the dashed line - in 

the soft + hard pomeron model) 
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For a quantitative description of the different features of the diffraction processes, 
we presumably need to take into account many different effects and a more com- 
plicated structure for the non-linear equation. 

We have shown that the most usual unitarization schemes could be recast into 
differential equations which are reminiscent of saturation equations I12j . Such 
an approach can be used to build new unitarization schemes and may also shed 
some light on the physical processes underlying the saturation regime. 
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